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Abstract. We define and study pseudo-diff'erential operators on a class of 
fractals that include the post-critically finite self-similar sets and Sierpinski 
carpets. Using the sub-Gaussian estimates of the heat operator we prove that 
our operators have kernels that decay and, in the constant coefficient case, are 
smooth off the diagonal. Our analysis can be extended to products of fractals. 
While our results are applicable to a larger class of metric measure spaces 
with Laplacian, we use them to study elliptic, hypoelliptic, and quasi-elliptic 
operators on p.c.f. fractals, answering a few open questions posed in a series 
of recent papers. We extend our class of operators to include the so called 
Hormander hypoelliptic operators and we initiate the study of wavefront sets 
and microlocal analysis on p.c.f. fractals. 

1. Introduction 

In this paper wc define and study pseudo-difFerential operators on metric mea- 
sure spaces endowed with a non-positive self-adjoint Laplacian such that the heat 
operator satisfies sub-Gaussian estimates. The motivating examples for our work 
are metric measure spaces constructed from p.c.f. fractals ([21. 38, 35]). However, 
our results capture the known facts from classical harmonic analysis on R" and 
Riemannian manifolds and can be used to define pseudo-differential operators on 
other metric measure spaces. For example, the operators defined in [16] fall under 
the hypothesis of our main theorems. 

Analysis on fractals from either the probabilistic or the analytic viewpoint has 
been the focus of intense study recently (see, for example, [15, 4, 21, 40] and refer- 
ences therein) . Several recent papers have studied properties of spectral operators 
on fractals. For example, [39] shows some new convergence properties of Fourier 
series on fractals with spectral gaps and establishes a Littlewood-Paley inequality 
for such fractals. In [20] a kernel formula for the resolvent of the Laplacian on any 
p.c.f. fractal is given; this generalizes Kigami's result for the Green function on 
p.c.f. fractals ([21]). These results have been extended to infinite blowups of frac- 
tals in [24] . Numerical results suggest that this result might hold for other spectral 
operators on fractals [1, 7]. The main result of [19] says that if a spectral operator 
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on a p.c.f. fractal is given by integration with respect to a kernel that is smooth 
and has a specific type of decay off the diagonal, then it is a Calderon-Zygmund 
operator in the sense of [31, Section 1.6.5]. In particular, the authors show that 
the Riesz and Bessel potentials and, more generally, the Laplace type operators are 
Calderon-Zygmund operators. On the other hand, in [11, 28], the authors show 
that spectral multiplier operators on metric measure spaces and products of such 
spaces are bounded on as long as the heat operator satisfies specific estimates. 
The Laplacian defined on some p.c.f. fractals and some highly symmetric Sierpinski 
gaskets satisfies these assumptions. 

We begin by presenting the standing hypothesis of our results in Section 2. We 
also briefly review the definition and main properties of p.c.f. fractals and the 
Laplacians defined on them. 

The pseudo-differential operators with constant coeflacients that are studied in 
Section 3 generalize the class of spectral multipliers studied in the papers mentioned 
above. We prove that these operators satisfy the symbolic calculus and that they 
are given by integration with respect to kernels that are smooth and decay off the 
diagonal, extending some of the results of [19]. Moreover, the class of pseudo- 
differential operators of order on p.c.f. fractals are Calderon-Zygmund operators 
and thus extend to bounded operators on L*, for all 1 < g < oo. In this context 
we therefore recover the results of [11, 28]. We extend our analysis to products of 
such spaces in Section 5. 

We define Sobolev spaces on spaces built out of p.c.f. fractals in Section 4 and 
prove that pseudo-differential operators with constant coefficients are bounded on 
them. We study elliptic and hypoelliptic operators on fractals in Section 6. Namely, 
we prove that a pseudo-differential operator satisfies the pseudo-local properties and 
that an elliptic operator is hypoelliptic. This gives positive answers to some open 
questions posed in [40, 6, 25]. An interesting class of operators that can be defined 
on fractals with spectral gaps are the so called quasielliptic operators ([6, 28]). 
We show that every quasielliptic operator is equal to an elliptic pseudo-differential 
operator, though there are quasielliptic differential operators which are not elliptic 
as differential operators. 

In Section 7 we extend the class of pseudo-differential operators to include op- 
erators for which the derivatives of the symbols have a slower rate of decay. As an 
application we show that the Hormander type hypoelliptic operators are hypoellip- 
tic. This extends one side of the classical result ([42, Chapter III, Theorem 2.1]). 
The converse is false in general, as is exemplified by the quasielliptic operators (see 
Subsection 6.1). 

Section 8 introduces the wavefront set and microlocal analysis on products of 
compact spaces built out of fractals. We show that pseudo-differential operators 
may decrease the wavefront set and that elliptic operators preserve the wavefront 
set, extending results from classical harmonic analysis (see, for example, [29]). We 
also describe the wavefront set for a few specific examples. 

In the last section, we study properties of pseudo-differential operators with 
variable coefficients on compact fractafolds. The main results of the sections are the 
continuity off the diagonal of the kernels of such operators and the L'-boundedness 
of these operators. These results cannot be obtained using the methods of [11] 
and [28]. We also conjecture that the pseudo-differential operators with variable 
coefficients of order are Calderon-Zygmund operators. 
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Some of the results proved here are extensions of the corresponding results from 
classical harmonic analysis (see, for example [42, 30, 32]) However the proofs of our 
results arc very different. The main reason for this difference is that the product 
of smooth functions is, in general, no longer in the domain of the Laplacian [5]. 
Therefore techniques that arc essential in real analysis like multiplication with a 
smooth bump arc not available to us. Wc frequently use the Borcl type theorem 
proved in [26] to decompose a smooth function in a sum of smooth functions. 

Acknowledgement. The authors would like to thank Camil Muscalu and Alexander 
Teplyaev for many useful conversations during the preparation of the manuscript. 
They also thank the referee for various suggestions, including pointing out the 
relevance of [16] to this topic. The authors express their thanks to Yin- Tat Lee for 
pointing out an error in an earlier draft of this article. 



2. Background 

For the main results of Sections 3 and 5 we need a metric space {X, R) with 
a Borel measure ^ and a negative self-adjoint Laplacian A. Wc assume that X 
satisfies the doubling condition, that is, there is a constant C > such that 

(2.1) KB{x, 2r)) < Cfi{B{x, r)) for all x e X and r > 0. 

We also assume that the heat operator e*'^ has a positive kernel ht{x,y) that 
satisfies the following sub-Gaussian upper estimate 

(2.2) htix,y)<c^t-^c^p(^^C2(^^^^^j^yy 

where ci,C2 > arc constants independent of t,x and y. In this expression both 
d and 7 arc constants that depend on X and (3 = d/{d + 1). Moreover we assume 
that hz{x,y) is an holomorphic function on {Rcz > 0}. 

Examples of spaces that satisfy the above hypothesis include the so called p.c.f. frac- 
tals and the highly symmetric Sierpinksi carpets. We review some of the definitions 
and properties of the p.c.f. fractals and of the Laplacians defined on them. For 
more details the reader can consult the books [21] and [40]. Recall that an iterated 
function system (i.f.s.) is a collection {Fi, . . . , i^jv} of contractions on W^. For such 
an i.f.s. there exists a unique self-similar set K satisfying (see [18]) 

k^Fi{k)\J---\Jfn{k). 

For oji, aj2, • • - (^n G {1, • • ■ ,N}, lu = luiuji . . .a;„ is a word of length n over the 
alphabet {1, . . . , N}. The subset K^, = Fu{K) F^^^ o ■ • • o F^^ (K) is called a 
cell of level n. The set of all finite words over {1, . . . , A^} is denoted by W^*. Each 
map Fi of the i.f.s. defining K has a unique fixed point Xi. We say that K is a 
post- critically finite (p.c.f.) self-similar set if there is a subset Vq C {xi, . . . ,xn} 
satisfying 



(K) fl F^, (K) C F^ {Vo) fl F^, (Fq) 

for any uj ^ lo' having the same length. The set Vq is called the boundary of K and 
the boundary of a cell is F^(Vb). One defines Vi = IJ- Fi{Vo), and, inductively, 
Vn = Ui Pi{yn-i) for n > 2. The fractal K is the closure of ij„ K- 
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The Laplacian on p.c.f. fractals may be built using Kigami's construction [21] 
from a self-similar Dirichlet energy form £ on K with weights {ri, . . . , rjv}: 

N 

1=1 

The existence of such forms is non-trivial, but on a large collection of examples they 
may be obtained from the approximating graphs as the appropriate renormalized 
limit of graph energies ([21, 40]). 

The second ingredient is the existence of a unique self-similar measure 

N 

where {/ii, . . . , /xjv} are weights such that Q < < 1 and ^ = 1, see [18]. 
Then the Laplacian is defined weakly: u G doniA^ with A^m — f ii 

£{u, w) = — / fvdji 
J X 

for all V G domf with v\vo = 0. The domain of the Laplacian depends on the 
assumptions that one makes about /. Kigami ([21]) assumes that / is continuous, 
but in this paper it will be more natural to assume that / is in i^(/i), which gives 
a Sobolev space (see Section 4). We write u G dom^p A if / G LP{^). 
The effective resistance metric R{x,y) on K is defined via 

R{x,y)^^ — mm{£{u) : u{x) = and u{x) = 1}. 

It is known that the resistance metric is topologically equivalent, but not metrically 
equivalent to the Euclidean metric ([21, 40]). 

The unit interval and the Sicrpinski gasket ([21, 40, 13, 5, 23, 36, 43]) are im- 
portant examples of p.c.f. fractals. Other examples are the affine nested fractals 
[12] that in turn are generalizations of the nested fractals [22]. 

Some of the spaces that we consider in this paper are built from p.c.f. fractals as 
in [35, 36]. In those papers the author defines fractal blow-ups of a p.c.f. fractal K 
and fractafolds based on K. The former generalizes the relationship between the 
unit interval and the real line to arbitrary p.c.f. self-similar sets, while the latter 
is the natural analogue of a manifold. Let w G {1, . • . ,N}°° be an infinite word. 
Then 

The fractal blow-up is 

GO 

X I I . ..F^^K. 

m— 1 

If C is an n cell in K, then F~^ . . . F~^C is called an (n — m) cell. The blow-up 
depends on the choice of the infinite word w. In general there are an uncountably 
infinite number of blow-ups which are not homeomorphic. In this paper we assume 
that the infinite blow-up X has no boundary. This happens unless all but a finite 
number of letters in w are the same. One can extend the definition of the energy £ 
and measure fi to X. The measure fj, will be a-finite rather than finite. Then one 
can define a Laplacian on X by the weak formulation. It is known that, for a large 
class of p.c.f. fractals, the Laplacian on an infinite blow-up without boundary has 
pure point spectrum ([43], [27]). 
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One can view a fractal blow-up as a collection of copies of K that are glued 
together at boundary points. Every point of X has a neighborhood which is home- 
omorphic to a neighborhood of a point in X . A fractafold [36] based on K is defined 
to be a set which satisfies this latter property. We will work on a restricted class of 
fractafolds, which may be thought of as more like triangulated manifolds. Specifi- 
cally we will consider fractafolds X that consist of a finite or infinite union of copies 
of K glued together at some of the boundary points. The fractafold X is compact if 
and only if we consider a finite number of copies of K. We suppose in the following 
that all the copies of K have the same size in X. If all the boundary points of the 
copies of K are paired, then the fractafold X has no boundary. When K is the 
unit interval this construction produces the unit circle. The next simplest example 
is the double cover of the Sicrpinski gasket, where one consider two copies of the 
fractal with corresponding boundary points paired. One can extend the definition 
of energy and Laplacian from K to a, fractafold based on K . An explicit description 
of the spectral resolution of the fractafold Laplacian for certain infinite fractafolds 
is given in [33]. 

Products of fractals provide another important class of examples for our results. 
An important point to keep in mind is that product of p.c.f. fractals is not a 
p.c.f. fractal. Strichartz described in [38] how one can extend the definition of the 
Laplacian and energy to products of fractals. 

The estimates (2.2) are known to be true for < i < 1 on a large number of 
p.c.f. fractals ([12], [4], [14], [15], [41]) and Sierpinski carpets ([2, 3]), with R being 
the resistance metric, the constant d being the Hausdorff dimension with respect 
to the resistance metric, and 7 a constant specific to the fractal. For the Dirichlet 
Laplacian on a p.c.f. fractal the estimate for t > 1 is immediate because the heat 
kernel decays exponentially at a rate determined by the smallest eigenvalue. In 
the case of the Neumann Laplacian one should instead modify the heat kernel by 
subtracting its projection onto the zero eigenspacc of constant functions, which has 
exponential decay controlled by the first non-zero eigenvalue. This latter change 
does not affect the definition of pseudo-differential operators because their action 
is on the eigenspaces with non-zero eigenvalues; accordingly we will abuse notation 
and refer to this modified heat kernel as just the heat kernel. No modification is 
needed in the case of blowups of p.c.f. fractals. For many fractals, lower bounds on 
the heat kernel are also known, but we will not need these in this paper. We also 
note that the heat kernel is holomorphic on {Rc z > 0} (sec [21]). 

We mention that estimates of the form (2.2) are also known for sums of even 
powers of vector fields on homogeneous groups [16]. However, we will not pursue 
these spaces in this paper. 

In this paper we write A{y) < B{y) if there is a constant C independent of y, 
but which might depend on the space X, such that A{y) < CB{y) for all y. We 
write A{y) ^ B{y) if A{y) < B{y) and B{y) < A{y). If f{x,y) is a function on 
Xi X X2, then we write Ai/ to denote the Laplacian of / with respect to the first 
variable and A2/ to denote the Laplacian of / with respect to the second variable; 
repeated subscripts indicate composition, for example A21 = A2 o Ai. We say that 
a function u is smooth if u € Dom A" for all n > 0. 
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3. Symbols and pseudo-differential operators on fractals 

In this section {X, R) is a metric measure space with a Laplacian A that satis- 
fies the conditions (2.1) and (2.2). We write P{X) for the spectral resolution of the 
positive self-adjoint operator —A. Examples of spaces for which the results of this 
section apply include compact or infinite blow-ups without boundary or products of 
copies of the same fractafold. In the case that X is either a compact fractafold with- 
out boundary or an infinite fractafold without boundary for which the Laplacian 
has pure point spectrum ([43, 27]), we write Px for the spectral projection corre- 
sponding to the eigenvalue A. If X is compact or A has pure point spectrum we fix 
an orthonormal basis { 0„ or { }ni^i. of L^{p) consisting of cigcnfunctions 

with compact support and write 2? for the dense set of finite linear combinations 
with respect to this orthonormal basis. If X is a product of such fractals then there 
is a natural basis of L^{fJ.) obtained by taking products of eigenfunctions on each 
fiber ([38]). The following are the main objects of study in this paper. 

Definition 3.1. For fixed m G K we define the symbol class S'™ to be the set of 

p G C°°((0, oo)) with the property that for any integer fc > there is Cfc > such 
that 



for all A > 0, where d is an in (2.2). 

Remark 3.2. The rationale for dividing m by d -I- 1 is that the Laplacian behaves 
like an operator of order d + 1. 

If p is any bounded Borcl function on (0, oo) then one can define an operator 
p(— A) via 



This operator extends to a bounded operator on (/i) by the spectral theorem. If 
p G S*™ with m > 0, then q{X) := (1 + A)~^p(A) is bounded and one can define 



Definition 3.3. For fixed m G M define the class '^DOm of pseudo-differential 
operators on X to be the collection of operators p(— A) with p G 5™. 

If X is a compact fractafold without boundary or A has pure point spectrum 
then the formula for a pseudo-differential operator is 



for p G 5™ and u ^ V, where A is the spectrum of —A. 

Proposition 3.4 (Symbolic Calculus). If pi G 5™^ and p2 G 5™^ then pip2 G 
^mi+m^ and 





p(-A) = (/-A) — g(-A). 



pi~A)u=Y,p{X)Pxu 



Pi(-A) op2(-A) =pip2(-A). 



Proof. Let fc > 1. Using Leibnitz's formula we have that 
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Therefore 



< ^ ('^)(1 + A)^(1 + A)^T 



Thus pip2 G 5'"!+'"^ Using the identity P{X){p{-A)u) = p{X)P{X){u) we have 
that 



Pi(-A)op2(-A)«= / pi(A)P(A)(p2(-A)ii) 

Pi(A)p2(A)F(A)H - (pip2)(-A)u. 



□ 

The main resuh of this section is Theorem 3.6, which says that if p G 5*° then 
p(— A) is given by integration with respect to a kernel that is smooth off the diagonal 
and satisfies specific decay estimates off the diagonal. When X is a fractafold based 
on a p.c.f. fractal we then obtain from [19, Theorem 1] that p{—A) is a Calderon- 
Zygmund operator in the sense of [31, Section 1.6.5]. We begin with a technical 
lemma. 

Lemma 3.5. Let a > and R > be fixed. Then 

(3-1) / cxp(-ci?(''+i)^2"^(l + C')-^)^^-i^dC< Ci?-" 

provided j > Ijgtll + 1. 
Proof. Consider the integral 



/ = / cxp(-^(l + e^ 
Jm 



2) exp(-A(l + a-"*')(^^rfC. 



Since 2(1 + > (1 + on (0, oo) we have 

for j > i. We use this bound for A > 1 and the obvious bound by C{j) for A < 1. 
Let A ^ cR'^'i+i)j2^n and no be such that A > 1 iff > uq. The series in the 
statement of the lemma is bounded by 

X ^ Tin (c' + l)(l-2j)'y X ^ f a (2j-l)7\ 

C{j) ^^+C{jn)R 2'H3+T-HT^j 

which converges for j as in the lemma. The estimate follows from the fact that 
2"o 2i □ 

c 

Theorem 3.6. Let p : (0, oo) C be an -symbol, that is p is smooth and for all 
k > there is Ck > such that 

Qk 



(3.2) 



<Cfc. 
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Then p{—A) has a kernel K(x,y) that is smooth off the diagonal of X x X and 
satisfies 



(3.3) 

and 

(3.4) 



\Kix,y)\<Rix,yr'' 



|A^A^/^(x,y)|<i?(x,y)-'^-('+^)(''+i). 



Proof. We begin with the Littlewood-Paley dyadic decomposition from [32. page 
242]. Let 7] be a C°° function with r?(A) = 1 if |A| < 1 and r?(A) = if |A| > 2 and 
let 5(A) = r/{X) - ?7(2A). Then supp(5 C {i < |A| < 2} and 

^<5(2-"A) = l, 

where for each A there are only two nonzero terms in the above sum. Then S is 
and we let Dk > for each fc > such that 



(3.5) 



d^ 



For n e Z let p„(A) = p(A)(5(2-"A). Then suppp„ C [2"-i,2"+i] (since we 
assume that p is defined on (0, oo)) and 



(3.6) 



5(A) = 5]p„(A). 



Moreover we can use (3.2) and the fact that the support of p„ has support in 
[2"~^, 2"+^] to bound -jrkPn by a constant that depends only on k: 



d'' 
dA^ 



dXi^-i 



3=0 



(3.7) 



<2-"^■C(fc)E|J-2"■'" = Cfc2-"^ 



2nj 



where C(fc) = maxC^Dj^k and Cfe = 2'=+iC(A:). 

Fix now 71 G Z and set /,i(A) = p„(2"A)e'^. Then supp/„ C [^,2] and for all 
fc > we have that 



dX''- 



d' 



/„(A)=5: J2"^— p„(2"A)e^ 



dX^' 



Using (3.7) and the fact that A < 2 wc obtain 



(3.8) 



d^ 



3=0 
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It follows immediately that there are constants Bk independent of n such that 
the Fourier transform fn{£,) of /„ satisfies \£,'^,fn{0\ ^ Bk, and thus for each fc G Z 

(3-9) lum < j^f^- 

for some constants Dk independent of n. 
Now by Fourier inversion 

so that 

(3.10) Pn{~A)u = /„(C)e^^'-^«^^«dC 
and the kernel of p„ (— A) is given by 

(3.11) K^{x,y) = / UOh^^M^,v)dt 

where hz{x,y) is the complex heat kernel, which is holomorphic on {Rez > 0}. A 
proof identical with that of Lemma 3.4.6 of [8] shows that there is C > such that 

(3.12) |/i,(x,j/)| <C(Rez)-^. 

Using the above estimate together with the heat kernel estimates (2.2), Lemma 9 
of [9] implies that 

(3.13) \h,{x,y)\<C2—{\z\cos{e))- — e 1^1 

where z = r(cos9 + ism9). For z = ^ — wc have that \z\ — 2~"(1 + ^^)5, 
cosfl = {l + C'^yi, and \z\ cos9 = 2"". Thus 

(3.14) \h^_,^^{x,y)\ < C2^2^exp(-^i?(a;,y)(''+i)''2"''(l+e')"^ 
and 

\Knix,y)\ < C{d,j)2^ J exp(-^i?(x,y)('^+i)^2"^(l + ^2)-^) (JTe^^^' 

for all J > 1, where the constants C{d,j) depend only on d and j. Lemma 3.5 with 
a = d and R = R{x,y) implies that if we choose j large enough then K{x,y) = 
J2n€Z ^ri{x,y) is defined and continuous off the diagonal and it satisfies (3.3). 
Finally, using (3.6), we see that K{x,y) is the kernel of p{—A). 

Next we want to prove the K{x, y) is smooth off the diagonal and show that the 
estimates (3.4) hold. For simplicity we give the complete argument only in the case 
I = and fc = 1, i.e. we show that 

(3.15) \AyK{x,y)\<R{x,y)-^'-\ 

Let g(A) = Ap(A) and (?n(A) = Ap„(A) for all ti G Z. By (3.7) we obtain that 
\qn\\)\ < C(fc)2(i-'=)" for all k > 0. If we set g„(A) 2^"g„(2"A)e^ one then 
obtains |3i'^''(A)| < A{k) where A{k) are constants independent of n. Then a 
computation similar to (3.9) shows that 

(3-16) \9nm < 



10 



M. lONESCU AND L. G. ROGERS AND R. S. STRICHARTZ 



for all A; > 0, with B(k) independent of n. Using the Fourier inversion formula we 
have that 

and a similar formula holds for (7„(— A). Hence 
and so 

\AyK,,{x,y)\ < C(d)2"2— J |5„(0|exp i-c{j)^0^^^\ d( 

^, , ,.. "(2^+1) f ( , ,i?(a;,y)(''+i''^2"''\ 1 
< C(,.,)2^ J cxp (^-»(„^-lL-^ j 

for all J > 1. Using Lemma 3.5 with a = 2(i + 1 wc find for j large enough that 
AyK(x,y) — AyKn{x,y) is well-defined and continuous off the diagonal and it 
satisfies (3.15). 

To prove (3.4) for general I and k, one can repeat the above steps for the functions 
g(A) = A'+V(A), g„(A) = A'+*^>„(A), 5„(A) = 2-"('+'^-)<7„(A)e^ and apply Lemma 
3.5 with a = d+ (^ + fc)(d+ 1). □ 

Corollary 3.7. Assume that X is a fractafold without boundary or a product of 
such fractaf olds. If p G S'^ thenp{—A) is a Calderon-Zygmund operator and, thus, 
it extends to a bounded operator on L'^ifJ.) for all 1 < q < oo and satisfies weak 
1 — 1 estimates. 

Proof. Using the estimates (3.3) and (3.4) (with I = and k = 1), Theorem 1.1 of 
[19] implies that p(— A) is a Calderon-Zygmund operator. □ 

Remark 3.8. The boundedness of p(— A) on i^(/i), 1 < q < oo, can be obtained 
also using the results of [11]. 

Corollary 3.9. If p £ S"^ then p(— A) is given by integration with respect to a 
kernel Kp that is .smooth off the diagonal. 

Proof. If p e S*™ then g(A) = p(A)(l + A)"™ G S° and Theorem 3.6 implies that 
q{—A) has a kernel Kq that is smooth off the diagonal. Then 

Kp{x,y) = iI-A^rKq{x,y) 

is smooth off the diagonal and it is the kernel of p(— A). □ 

As a consequence, we obtain that, for Res > 0, the Bessel potentials 

(3.17) (/- A)-^w= ;^(1 + A)-^Pau, 

AeA 

and the Riesz potentials 

{~A)-'u = A"-' -Paw, 
AeA 

have smooth kernels and are bounded on LP{^), for I < p < oo, with operator norm 
of at most polynomial growth in Ims (when Res = 0). These facts were previously 
proved in [19]. 
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4. SoBOLEv Spaces 

In this section wc assume that X is a compact fractafold without boundary based 
on a p.c.f. fractal K ^ an infinite blow-up of K without boundary, or a product of 
copies of the same fractafold without boundary. Moreover, we assume that if X 
is non-compact then A has pure point spectrum. Therefore we can consider that 
L'^in) is spanned by compactly supported eigenfunctions of —A. As before, we fix 
an orthonormal basis {(/>„ } of L'^{^j) consisting of such eigenfunctions and we write 
2? for the space of finite linear combinations of 0„'s. 

Definition 4.1 (L^-Sobolev spaces). For s > and w £ define 

Wnfu. ■.^Y.^l + \)^^\\PM\l 
AeA 

We say that u G iJ* if and only if ||w||^fs < oo. 

For s > if M G then u G L^(/i). So for s > there is no harm in starting 
with u 6 in Definition 4.1. Moreover if si < S2 then H"^ C H'^^ . 

For s < we define using the distribution theory developed in [25] . Namely, 
for s > 0, H^'^ is the dual of H'^ via 

{f,f) = ^{.f,4>n){f,(f>7i) 

for / e H^'^ and ip E H\ The fact that the above linear functional is bounded is 
a consequence of the Cauchy-Schwartz inequality: 

n 

Lemma 4.2. For s a positive integer, u G if and only if{I — A)''u G L^{n) 

for all k < s. 

Proof. Notice that 

AeA 

From this equality the statement follows. □ 

Proposition 4.3. If p G S"" then p{-A) : H'-'"- for all s and m. 

Proof Recall that Px{p{-A)u) =p{\)Pxu. Then 

|b(-A)^.||?,._ = E(l + A)W||P,(K-A)u)||2 
AeA 

^J2(^ + xf-^\pi\)\'\\Pxu\\l 
AeA 

<Y.{i+x)^\\Pxu\\i = \\urH.. 

AeA 

Thus ^(-A)^ G ff"-'" whenever uE H". □ 
Proposition 4.4. For s > 0, H'^ equals the image of L^{^) under (/ — A)^''/('^+^). 
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Proof. Let uje L^{fi) such that u = {I - A)-''/('^+i)/. Then 



Therefore 



~ (1 + \y/(d+i) 



AsA ^ ' 

Now let ueH'' and define f = {I - Ay^^'^+^^u, that is 
Then 

AeA 

So / e L2(/i) and u = (/- A)-''/('^+i)/. □ 

We extend next the definition of Sobolev spaces to spaces, using Proposition 
4.4 as the starting point. Namely, we are going to replace L'^in) with LP{^) (see 
[30]). 

Definition 4.5. We define the Sobolev spaces Lf for s > and f < p < oo to 
be the images of under (/ — A)~^/(''+^^ with the norm 



||(/-A)-«/(^'+i)/|| 



Note that by Corollary 3.7 we may regard Lf as a closed subspace of L^, with 
< c||u||^p. For s = we have Lq = and for p = 2 we have = H\ 

Proposition 4.6. // < sq < si < oo and 1 < po,pi < oo then the complex 
interpolation space [-^fo'j-^fjjg may be identified with if where < 9 < 1, s = 
(l-0)so + esi, and ^ = ^ + 

Proof. The proof follows as in the Euclidean case using the comments following 
Proposition 3.7. □ 

Proposition 4.7. Ifp G S"™ thenp{~A) is a bounded operator from if into L^^„^. 

Proof Let / e Lf. Then there is g G LP such that f ^ {I - A)-'/'^'^+^^ g. Then 
p(-A)/ = p(-A)(/ - A)-*/(''+i)g. By Proposition 3.4 and Corollary 3.7 

(/ - A)('*-"')/(''+i)p(-A)(/ - A)--^/(''+i',9 e pp. 

This is equivalent with p{-A)f = (I - A)-("-'"'/(''+i)g e Ll_^. □ 

Lemma 4.8. The Sobolev space L^^^j^^ equals domLp{A) . 

Proof. Let u E domLp{A). Then u G dom{£) and there is / G such that 
/ = -Alt. That is 

fvdjjL = £{u,v) for all v G dom(£). 

li V = 0Aj where (px is any eigenfunction corresponding to A, then A ^ u(j)\d^ = 
/ /^Acf/i- This implies that (/ - A)"i(/ + u) = u. 
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For the converse, notice first that if Ai,A2 G A and (px- is any eigenfunction 
corresponding to Ai, i = 1,2, then £{4>\^,(j)\^) = Ai(S(Ai, A2). Let u = [I — ^)^^g 
for some g G L^. Set f = g — u E . For any A e A wc have 

f(t)\d^i, = j {g - u)(t)xdfi = J g(t)xdfi. 

Also 

It follows by linearity and density that J fvd^ ~ £{u,v) for all v e dom£. Thus 
-Am = /. □ 

Theorem 4.9 (Sobolev embedding theorem). If s < ^ then Lf C for | = 

Proof. We have all the ingredients needed to use the same proof as in [37, Theorem 
3.11]. □ 

5. Pseudo-differential operators on Products of Fractals 

In this section we extend the definition of pseudo-differential operators to prod- 
ucts of metric measure spaces that satisfy the hypothesis of Section 2. Let N > 2 
be fixed and let (ATi, . . . , {Xn,Rn) be N metric spaces such that Xj has mea- 
sure and Laplacian A^. In case that Xi is a fractafold, then we consider that 
Ai is defined using a self-similar Dirichlet energy £i. We assume that the heat 
kernel /i'-*' associated to A^ satisfies the estimates (2.2) with d ^ di, j = ji, 
R = Ri, for alH = 1, . . . , A^. Let X = Xi x ■ • ■ x Xn be the product space and let 
fj, = X ■ ■ ■ X /ijv be the product measure on X . We write x = {xi, . . . , Xn) for 
elements in AT, A = (Ai, . . . , Xn) for elements in (0, cxd)^ and A = (Ai, . . . , A^v). 
Recall that there is a unique spectral decomposition Px such that 

(5.1) PxUi X U2 X ■ ■ ■ X UN = Pxi PX2 ■ ■ ■ ® P\n'^1 X • • • X UN, 

where mi x • • • x Un{x) ~ Ui(xi)u2{x2) ■ ■ ■ un{xn) and 

Pxt(E)Px2 ■ ■ ■®-PA„wi(a;i)M2(a;2) • ..unIxn) = PxiUi{xi)Px2U2{x2) ■ ■ ■ Pa„ww(2^7v), 

Pxi is the spectral projection corresponding to the eigenvalue Ai of — A^. 

Recall from [38] that if {<pn} is an orthonormal basis for i^(/Ui) consisting of 
compactly supported eigenfunctions of — Ai, then 

<Afci,...,fe«,(a^) = 0fei(^i) <^f„(a^w)> 
ki € Z for alH = 1, . . . , A^, form an orthonormal basis for i^(/i). In this case we 
write 2?^ be the set of finite linear combinations of {4>ki....,kN}- 

Definition 5.1. For m G Z we define the symbol class 5™ on X to consist of the 
set of smooth functions p : (0, 00)^ — > C that satisfy the property that for every 
a = («!, ■ • ■ , ckAr) £ there is a positive constant Ca such that 



(5.2) A" 



<C„(1 + A)^, 



where A" = X'^'X^^ • ■ • A^", and 



^p(A) = ^...^p(Ai,...,A^). 
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Remark 5.2. Note that each Xi is positive and we omit the absolute values of A 
in the above definition. Notice also that if p satisfies the condition (5.2) then it 
satisfies also the condition 



(5.3) 



|A| 



|q| 



dX' 



< C„(l + A) — 



where |A| — \/ Xf + ■ ■ ■ + Xj^ if A e (0, oo)^, and \a\ = ai + • ■ • + un, if a G 
N^. Condition (5.2) is usually called the Marcinkiewicz condition while (5.3) is 
the Hormander condition. Many authors define psuedo-differential operators on 
Euclidean spaces or manifolds using the Hormander condition. We chose to use 
Marcinkiewicz condition in the definition because it is more general and it makes 
the proof of the main theorem in this section more transparent. 

If p : (0, oo)^ — > C is a bounded Borel function then we can define an operator 
p(— A) acting on L^{fJ.) via 



(5.4) 



P{-A)U : 



pWPxiu), 



(0,oo) 



where Pa is defined via (5.1). 

Definition 5.3. For m £ Z define the class ^DO^ of pseudo-differential operators 
on X to be the collection of operators p{—A) with p G S*™. 

The spectral theorem implies that, if m = 0, then p{—A) extends to i^(/i). 
Moreover we will show that, as in the single variable case, these operators are given 
by integration with respect to kernels that are smooth off the diagonal and, if X is a 
product of fractafolds, they are Calderon-Zygmund operators on X. The following 
is the main theorem of this section. 



Theorem 5.4. Suppose that p G S'°, that 



is, p is smooth an 



(5.5) 



A" 



for all a = (ai, • • • , un) <= N^. Then p{— A) , where A = (Ai, . . . , Ajv), is given by 
integration with respect to a kernel Kp that is smooth off the diagonal and satisfies 

N 

(5.6) \Kp{x,y)\<l[Rkixk,ykr'''' 



k=l 



and 
(5.7) 



N 



A%Al^^^Kp{x,y) < i?,(x,2;,)-ft(*+i)i?,(:r„y,)-''^^'^+'^ H Ri^k^yk)-"" , 

k=l 

for all i,i = 1, . . . , A^, where x = (.ti, . . . , xat), y = (j/i, . . . , y^r) and A^^i is the 
Laplacian with respect to Xi and Ay j- is the Laplacian with respect to yj. 

Proof. For simplicity we prove the result for N ^ 2. The difference between this and 
the general case is a matter of notation. We use the Littlewood-Paley decomposition 
from [32, page 242] in each variable. Let S be the smooth function from the proof 
of Theorem 3.6 and let Dk be the bounds from (3.5). For n,m € Z define 

P„,™(A) ==p(A)<5(2-"Ai)5(2-'"A2), 
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where A ~ (Ai, A2) G (0, 00)^. Then pn.m is a smooth function with 
(5.8) suppp„,™ C [2"-\ 2"+i] X [2"-\ 2™+i], 

for aU n, m G Z. 

Let n, m G Z be fixed. We use Leibnitz formula to estimate the bounds on the 
derivatives oi Pn,m- If a = {a\^a.2 

) G W- then we have 
^Pn,™(A) = ^^(p(A)5(2-"Ai)<5(2-"A2)) 



a"i(5(2-"A 



j=0 k=0 



Q!2 

A: 



i 



a2 
k 



(9V(A) 

^aT 



2-m(o!2-fc)^(a2-fc)(|2~"'A2) 



• (5("i-J)(2-"Ai),5("^-'=)(2-™A2) 
Therefore, using (5.5) we obtain that 



p(A) ) 2^"("i^^^2^™'"^~'^^ 



9" 

^Pn.m(A) 



"1 



a2 ^ C(,-fc)2"J2""-- 
AiA2'= 



< C(a)2-""i2-™"^^ 



3.k 



Oil 

i 



a2 \ 2"J2"''2^+'' 
k J 2^3 2™-^ 



ai 

i 



"2 \ 2^+'= 
A: 



(5.9) 



< C(a)2-""i2-"'°^^ 



where C(a) = maxC(j fc) and we used in going from the first hue to the second that 
Ai > 2"-i and A2 > 2'"-i. Define now /„,,„(A) = p„,™(2" Ai, 2™A2)e^i+^^ Then 



(5.10) 



supp/„,„(A) C 



If a = (ai,a2) G then using Leibnitz's rule one more time we obtain that 



^/".™(^) = £^(P«.™(2"Ai,2"A2)e^^+^^) 



9" 



9A 



"2 ^ ryrnk ( 9^ 
9A§ 



p„,™(2"Ai,2™A2) e^i+^- 



3 



0L2 

k 



j=0 k=0 

Thus, using (5.9) and (5.10) we obtain that 



2"^'2"'=^J-^^p„,™(2"Ai,2'"A2)e^^+^= 



^^/«,7ri(A) 



< 



(5.11) 



1Q1+Q2„4 



I ( 2"-'2"'=C(j.fe)2""-'2""'=e'' 



< C(a)2"i+"2e'' ==: 
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Therefore the Fourier transform of fn,m satisfies 



d' 



for all a S N^, where Ba > 0, and, thus, 



(5.12) 



\fn.,miO\ < 



Dr. 



for all ^ — (Cii ^2) G and a G N^, where are positive constants depending 
only on a. Using the inverse Fourier transform we have 

p„,™(2"Ai,2™A2)e^^+^^ = ^ J fnMOe'^'^d^A 

Therefore 

Thus by spectral theory we have 
Hence the kernel of p„.„i(— A) is 



for all X — {xi,X2),y = (1/1,2/2) & X, where ft-^'^ is the heat kernel corresponding to 
Ai, i = 1,2. Using inequalities (3.14) and (5.12) we obtain that 



(2) 



\K^,n^{x,y)\ < C(d,j)2<'i+i / cxp 



C7i?i(x,y)(''i+i)'^2"'^i' 



f ( C72i?2(x,j/)(''^ + l)T^2'"^^ ^ l_ 



■ 2<*2 + l 



exp 



2(1+^1)- 



(1+^1)'^- 



for all j, fc > 1. Lemma (3.5) implies that (if we choose j and k large enough) 



\Kp{x,y)\ 



<C(d)i?i(xi,yi)-'*ii?2(x2,y2) 



\-d2 



Repeating the above proof for the function '7/3i.^2(A) = Af^A^^p(A), we obtain that 
the kernel Kp{x,y) is smooth off the diagonal and satisfies the following estimates 



Af\A^^ifp(x,2/) <C(d,/3i,/32)i?i(x,y)-^^('^^+i)-'^^i?2(x2,y2)-^^('^^+i)-^^ 



□ 



Corollary 5.5. Assume that X is a product of fractafolds without boundary. If 
p € S'^ then p{—A) is a Calderon-Zygmund operator and, thus, it extends to a 
bounded operator on L'^{fJ.) for all 1 < q < 00 and satisfies weak 1 — 1 estimates. 

Proof. Our theorem 5.4 and Theorem 6.1 of [19] imply that p(— A) is a Calderon- 
Zygmund operator. □ 
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Remark 5.6. The boundedness of p(— A) on for all 1 < g < cx) can be deduced 

also from the results of [28]. 

Corollary 5.7. If p G S"^ then p(— A) is given by integration with respect to a 
kernel that is smooth off the diagonal. 

Example 5.8. Consider the Riesz potentials ([38, page 591]) 



Clearly Pi{Xi, . . . , Ajv) = aTPh^TJv ' * ^ 1, . . . are in 5° and thus the Riesz po- 
tentials arc bounded on LP{fj,), for 1 < p < cxd, and arc given by integration with 
respect to smooth kernels. 



In this section we specialize the set-up of section 5 to the case when X is the 
product of N fractafolds without boundary. If Xi is not compact, we assume 
that Ai has pure point spectrum. We still write A = {Xi, . . . , Xn) G (0,oo)^, 
X ~ {xi, . . . ,xn) G X and A = (Ai, . . . , Aat). We begin by extending the definition 
of elliptic operators from [25, Definition 8.2] to elliptic pseudo-differential operators. 
We show that any pseudo-differential operator satisfies the pseudo-local property 
and the converse inclusion is also true for elliptic operators. In other words, we prove 
that elliptic operators on fractals are hypoclliptic. These results give a positive 
answer to some open questions posed in [37, 25]. 

We use the theory of distributions on fractals as developed in [25] for our defi- 
nitions and results. Most of the results, however, are true for more general metric 
spaces that we consider in Section 2 if one replaces the word "distribution" with 
"function" in what follows. 

Recall that the space of test functions T>{X) consists of all smooth functions 
with compact support. The space of distributions on X is the dual space V^X) of 
2?(X) with the weak-star topology ([25, Definition 4.1]). A distribution u is smooth 
on the open set C X ([25, Definition 8.1]) if there exists a smooth function / on 

such that 



for all (j) G 2?(r2). If flu is the maximal open set on which u is smooth, then the 
singular support of u is ([25, Definition 8.2]) 

sing supp(m) ^ X \ nu- 

If p G S"" then one can extend the definition of p(— A) to the space of distributions 
via 



for aU (f> e V{X). 

Definition 6.1. We say that p(— A) satisfies the pseudo-local property if 




Px, (E) ■ ■ ■ (E) Px^u, i = l,...,N. 



6. Elliptic and Hypoelliptic operators 



(6.1) 




(p(-A)w,0) = (?i,p(-A)0) 



sing supp(p(— A)ii) C sing supp(u). 



Theorem 6.2. If p is in S*™ then p{~A) satisfies the pseudo-local property. 
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Proof. Let u e I)'{X) and let Ku = sing supp(u). Then u\a is smooth, where 
H, ~ X \ Ku. We need to show that Pu\fi is also smooth. Let K be be any finite 
union of cells in il and let il' the interior of K. We identify u\a with the smooth 
functions given by (6.1). Then, using [25, Theorem 7.7], we can find a smooth 
function ui such that ui\k = u\k, and ui|f2<= ~ 0. We write then u — ui + U2, 
where U2 is supported on K''. Then Pu = Pui + Pu2. Since ui is smooth it follows 
that Pui is smooth. Therefore we only need to prove that Pu2\iv is smooth. For 
this we arc going to use the structure theorem for distributions [25, Theorem 5.10] 
(see also the comments following Assumption 7.2 in the same paper). Then U2 may 
be written as a locally finite sum of the form U2 = J2 ^"^^Vj ot U2 = '^2 ^^"'^'^^fji 
where the i^j are Radon measures and the fj are continuous functions with compact 
support. Notice that the support of each i^j is a subset of K'^ and, respectively, the 
support of each fj is a subset of K'^. 

Assume first that U2 ~ J2 ^"^^ i^j where I'j arc Radon measures. Let (j> G T>{fl'). 
Then 



where Kp is the kernel of p{—A) provided by Theorem 5.4. The above expres- 
sion makes sense because the support of is a subset of fi', while the supports 
of i/j are subsets of A'^. Since Kp is smooth it follows that the function x t— !> 
J2 J A"^^ Kp{x,y)dvj{y) is smooth for all x S f2'. Thus p(— A)u2|o' is smooth. A 
similar proof shows that p(— A)?i2|o' is smooth if U2 = It follows that 

p(— A)w is smooth on 17'. Since K was arbitrarily chosen, we have that p{—/S.)u is 
smooth on Q,. □ 

We turn now to the study of elliptic pseudo-differential operators. The first part 
of the following definition is from [25, Definition 8.4]. 

Definition 6.3. An elliptic polynomial of degree m on is any polynomial q of 
degree m satisfying the property that there exist two positive constants c and A 
such that 



An elliptic differential operator on X is an operator (/(—A) for some elliptic poly- 
nomial (j(A). More generally, an elliptic pseudo-differential operator of order m on 
X is a pseudo-differential operator p(— A) whose symbol p G 5*™ satisfies (6.2). 

Remark 6.4. Notice that we require that condition (6.2) holds only for nonnegative 
values of Ai, . . . , Ajv. In particular, the polynomial p{\) = Ai + • • • + Ajv is elliptic 
according to our definition. 

Theorem 6.5. An elliptic pseudo-differential operator p(— A) is hypoelliptic, mea- 





3 




(6.2) 



|g(A)| > c|A|^ for all A G M+ such that |A| > A. 



ning that 



sing supp(p(— A)u) = sing supp(u). 
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Proof. Since p is elliptic the set of its zeros is compact and p{X) ^ for all A G 
with |A| > A. Let ^ G C°° be a cut-off function which is on a neighborhood of 
the zeros of p and is 1 for |A| > A. Define pi(A) = ^(A)/p(A), so that pi(A) is zero 
on a neighborhood of the zeros of p. Then p{X)pi{X) = 1 if |A| > A. Moreover 
pi G S~"\ since for |A| > A we have that (see. for example, [42. Proof of Theorem 
1-3]) 

Q^+Q^H ha^— CK 

Using Lemma 3.4 and Theorem 5.4 we have that 

Pi(-A)p(-A)u + Ru, 

where R G 5^°° has a smooth compactly supported symbol r(A). Hence R is 
infinitely smoothing. It follows that if il is an open set such that p(— A)u|q is 
smooth then u|si A)p(— A)u|n — Ru\n is smooth. Therefore sing supp(u) C 

sing supp(p(— A)u). Since the converse is provided by Theorem 6.2, it follows that 
p(— A) is hypoelliptic. □ 

6.1. Quasielliptic operators. It is known ([6, 10]) that for some fractals the set 
of ratios of eigenvalues of the Laplacian has gaps. That is, there are a < (3 such 
that A/A' ^ {oi, /3) for all A 7^ A' in the spectrum of —A. Consider the product 
X = Xi X X2 of two copies of the same fractafold without boundary that have a 
spectral gap Let Ai and A2 be the Laplacian on the fractafold viewed as 

acting on Xi and, respectively, on X2. As before, we write A ~ (Ai, A2). For a 
real number a in (a, /?) and e > define the cone 

ra,s = {(Ai,A2) gK^ : |Ai -aA2| <eA2}. 

If a G (a, /3) then there is e > such that Fq.e does not contain any pair of 
eigenvalues of — A. Let 25 be the collection of such cones. Denote by G the collection 
of cones in that do not intersect any cone in &. Thus any pair of eigenvalues of 
—A lies in some cone in Q. Moreover, any pseudo-differential operator will "live" 
on Q, as we will make it precise in Lemma 6.9. 

Definition 6.6. We say that an m-symbol p is quasielliptic if it satisfies (6.2) for 
all A G F and F G ^, with the constant A depending on F. A pseudo-differential 
operator p(— A) is called a quasielliptic operator if p is a quasielliptic symbol. 

Example 6.7. If a G (a,/?) then the operator Ai — aA2 is a quasielliptic operator 
([6; 28]). An intriguing consequence of the existence of spectral gaps is that, if 
a G (a,/?) then the operator l/(Ai — aA2) is bounded on LP{ii) for all 1 < p < cxo 
[28]. The next proposition show that the quasielliptic operators are, in fact, elliptic 
operators. 

Proposition 6.8. A quasielliptic pseudo-differential operator p{— A) is equal to an 
elliptic pseudo-differential operator. 

This proposition follows immediately from the following lemma which says that 
a quasielliptic symbol equals an elliptic symbol outside a cone contained in the 
spectral gaps. Note that even if p(— A) is a differential operator, the elliptic operator 
is only a pseudo-differential operator. 

Lemma 6.9. Let p be a quasielliptic m-symbol. Then there is an elliptic m-symbol 
p such that p{X) = p{X) for all X GT and allT G G- 
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Proof. Let a G (a, (3) and e > such that the cone Ta^s is a maximal cone in © 
and let F = \ Ta,e- Then F is a union of two simply connected cones. Since p 
is quasicUiptic, there are c > and A > such that |p(A)| > c|A| for all A S F with 
|A| > A. Since F^ := {A S F : |A| > yl} is a union of two disjoint simply connected 
sets, we can define the function q{X) = logp(A). Let 7'(A) and s(A) be the real part 
and the imaginary part, respectively, of q{X). By hypothesis, r(A) > -£p^ log |A| for 
all A € Ta- Let f(A) and s(A) be two smooth extensions of r(A) and s(A) such 
that f (A) > ^ log |A| for aU A G with |A| > A. Let p(A) = exp(f(A) + iS{X)) 
and extend it so that p(A) = p{X) for all A G F. Then p is the desired elliptic 
m-symbol. □ 

Proof of Proposition 6.8. Let q be an elliptic symbol q such that p{X) = q{X) for 
all A G F and all T G Q. Then q{~A) is an elliptic pseudo-differential operator. 
Moreover, since p{X) = q{X) for all A G Ai x A2 (where A,; is the spectrum of — A^), 
it follows from the definition that p{—A) = q{—A). Therefore p{—A) is equal to 
an elliptic pseudo-differential operator. □ 

As an immediate consequence of the above Proposition and Theorem 6.5 we 
obtain the following corollary, which answers an open question posed in [25]. 

Corollary 6.10. A quasielliptic operator is hypoelliptic. 

7. p-TYPE SYMBOLS AND PSEUDODIFFERENTIAL OPERATORS 

In this section we define and study pseudo-differential operators for which the 
derivatives of the symbols have a slower rate of decay. These operators still have 
a kernel that is smooth off the diagonal, even though they might not be bounded 
on i'^-spaces, as the well known examples in the Euclidian setting show. As an 
application, we study the so called Hormander type hypoelliptic operators and 
prove that they are hypoelliptic in this more general setting. In the following we 
use the set-up of Section 3. The definitions and main theorem of this section can 
be extended without difhculty to the product setting so we omit the details (see 
Remark 7.5). We will use the fact that the equivalent of Theorem 7.3 holds for 
products of measure metric spaces in Subsection 7.1. 

Definition 7.1. For fixed m G M and < p < 1 define the symbol class S™ to be 
the set of p G C°°((0, 00)) with the property that for any fc > there is Ck{p) > 
such that 



<Ckip){l + X)^ 

UA 

for aU A > 0. 

Definition 7.2. Define the operator class '^DO^ by 

/•oo 

p(-A)it= / p{X)P{X)u 
Jq 

for pe S"^ and uGV. 

Theorem 7.3. Let < p < 1 and let p : (0, 00) <C be an Sp-symbol, that is, 
for each fc > there is Ck{p) > such that 

Qk 



(7.1) 



< Cu{p). 
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Then p(— A) has a kernel K{x, y) that is smooth off the diagonal of X x X . 



Proof. The proof of this theorem follows in large the proof of Theorem 3.6. We 
will present the main steps of the proof skipping the computations that arc very 
similar. Let 5 be the function from the Littlewood-Palcy decomposition and let 
p„(A) =p(A)5(2-"A). Then suppp„ C [2"-i,2"+i] and 

'^'^ < 2-""'= for all fc > 0. 



P«(A) 



Define p„(A) = p„(2"''A). Then suppp„ C [2"(i-p)-i, 2"(i-'')+i] and for each fc > 
there is a constant Ck{p) > such that 



Let /„(A) =p„(A)e2"(i-'') . Since e^-d-"' < for all A e [2"(i-p)-i, 2"(i-'')+i]^ it 
follows that for each A; > there is a constant Ak{'y) > independent on n such 
that 

d'^fniX) 



< Ml)- 



Therefore 

le'^/JOl < 1 for allfc>0. 

In particular we have that for each fc > there is a constant Dk{'y) independent on 
n such that 



(7.2) 



|/„(2"(''~i)OI < .7^ ..w, for ah fc > 0. 



(l+^2)fe 

Using the Fourier inversion formula we have that 
p„(2"''A)e^ = ^ 

Therefore 



PniX) = ^ / fniOe- 



1 

2^ 



*^)d^ 

^^r^)d^, 



which, using the substitution u = ^/2"'^'' becomes 

2"{p-i) 



PniX) 



2n 



Thus the kernel of p„(— A) equals 

2"(p-i) 



Knix,y) 



27r 



Using the estimates (3.14) and (7.2) wc obtain that 

\Kn{x, y)\ < 2TlfT(rf+(l-P)('i+l)(2j-l)) 

.(-^i?(x,y)('^+ih2"''(i+a-"*') (^^d^, 



cxp I 
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for all j > 0. Lemma 3.5 implies that if we pick j such that 

{d+{l-p){d+l){2j-l){j+l) , 1 



j > 



then K{x,y) = '^n-^'ri{x,y) converges. We can pick a positive j in the equation 
above provided that p > ^j^^. Then it follows that the kernel K satisfies the 
estimates 

\K{x,y)\ < j/)"''<^+i)-i . 

A similar argument for q{X) = A^p(A), j > 1, shows that K is smooth off the 
diagonal. □ 

Corollary 7.4. // < p < 1 and p £ 5™ then the kernel of p{—A) is smooth 
off the diagonal. 

Proof. The proof is very similar with that of Corollary (3.9). □ 

Remark 7.5. Note that for p <1 the operators p(— A) are not Calderon-Zygmund 
operators and they might not extend to for all values of q. Also, one can easily 
extend the definition to product of fractals and adapt the proof of Theorem 5.4 to 
show that the kernel of a pseudo-differential operator of type p with ^^^-j- < p < 1 
is smooth off the diagonal. 

7.1. Hormander type hypoelliptic operators. Let X be a product of N mea- 
sure metric space as in the setting of Section 5. Recall that we write A = (Ai, . . . , An). 

Definition 7.6. ([17]) We say that a smooth map p : (0, cx))^ — > C is a Hormander 
type hypoelliptic symbol if there are e > and A > such that 



(7.3) 



< CalAp'^l"! for |A| > A, 



P(A) 

where Cq are postive constants for all a e N^. 

Theorem 7.7. Suppose that p is a Hormander type hypoelliptic symbol with e > 
^ijq-j-. Then p{— Is) is hypoelliptic. 

Proof. We know from (7.3) that the set of zeros of p is compact. Let be a smooth 
function that is on a neighborhood of the zeros of p and (/5(A) = 1 if |A| > A. Set 

q{\) = ip{X)p{X)-\ Then 

^'^w= E aA^^w-^^w-^w^-^ 

oiH ha^i 

if |A| > A. It follows that q G S°. Since q{X)p{X) = 1 if |A| > A we have that 

g(-A)p(-A) = u + Ru, 

where R is an infinitely smoothing operator. Theorem 7.3 implies that p(— A) is 
hypoelliptic. □ 
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8. Wavefront set and microlocal analysis 

In this section wc initiate the study of the microlocal analysis on fractals. Namely, 
we define the wavefront sets on fractals and provide a few concrete examples. As in 
classical harmonic analysis, the wavefront sets contain information about the sin- 
gularities of distributions. We show that a pseudo-differential operator reduces, in 
general, the wavefront set. If the operator is elliptic, then the wavefront set remains 
unchanged. This is just a beginning of this study. We expect that the wavefront set 
will play an important role in the future, but we can not expect results analogous 
to classical analysis on propagation of singularities for space-time equations being 
controlled by the wavefront set, as there are localized eigenfunctions which prevent 
singularities from propagating at all. We continue to assume the set-up and nota- 
tions of Section 6. Namely, X is the product of N fractafolds without boundary, /i 
is the product measure on X, and A = (Ai, . . . , A^r). In addition, we assume in 
this section that each fractafold Xi is compact. 

Definition 8.1. Let F denote an open cone in MJ^ and VL an open set in X. We 
use Lpa^ to denote normalized eigenfunctions corresponding to eigenvalues Aq^. , 
and set Aq ~ A^^ + ■ • • + Aa„. A distribution u is defined to be C°° in fi x F if 
it can be written on O as a linear combination of eigenfunctions with coefficients 
having faster than polynomial decay over the eigenvalues in F. More precisely, if 
there is a sequence 6„ and a function v with v\n ^ u that has the form 

(8.1) t' = Cgipa^ -S) ipa^ -S) ■ ■ ■ (E> ipaM i 

a 

for values Cq. such that |cq| < 5„(1+Aq)~"/^''+^-' for all n and all {Aqj , . . . , Aq,„ } G F. 
We define the wavefront set of u, WF(u), to be the complement of the union of all 
sets where u is C°°. 

Remark 8.2. If u is a smooth function on X then WF(m) is empty. More generally, 
sing suppu is the projection of WF(u) onto X. 

Remark 8.3. If the fractals Xi have gaps in the set of ratios of eigenvalues, then 
there are special cones F for which every u is C°° on AT x F because F contains no 
eigenvalues (sec Subsection 6.1). 

Proposition 8.4. If p e 5"" then MVFp{-A)u C WF(w). // in addition p(-A) is 
elliptic then WF(p(-A)it) = WF('u). 

Proof. Let m be a distribution on X and suppose that is an open subset of X 
and F is an open cone in such that u is C°° in fl x T. Let v be as in (8.1). 
Then Pv\n = Pu and 

Pv = ^ Cap{Xai , ■ • ■ , Aa„ )</?qi ® • • • «) , 

a 

with |c„p(Aa, , . . . , Aa„)| < 6„(l+A„)-("-™)/('^+i) for aU {A„, , . . . , A,„ } e F. Thus 
Pu is C°° in ri X F and, since and F were arbitrary, WF(Pu) C WF(u). 

If P is elliptic then wc can write, as in the proof of Theorem 6.5, u ~ PiPu + 
Ru, where Pi is a pseudo-differential operator of order — m and i? is a smoothing 
operator. The conclusion follows immediately from this. □ 
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Example 8.5. Let N ~ 2 and let u = mi (8)M2 be a tensor product of distributions. 
From the definition u is C°° on ((sing supp ui^ x (sing supp 1*2)'^) xMj, thus WF(u) 
is a subset of the complement of this set. It is also easy to check that 

(sing supp(ui) X sing supp(u2)) x 1^2^ C WF(u). 

Let xi S Xi and X2 G X2 and suppose that ui is smooth near xi but X2 G 
sing supp(u2). Then there is a smooth function vi on Xi and a neighborhood Qi 
of Xi such that ui = vi on f2i. Moreover we can write vi = Cjipj with the 
Fourier coefficients Cj having faster decay then any polynomial ([25, Theorem 3.5]). 
Therefore for any n £ N there is bn.j > such that 

|c,| < V,(l + A,)-". 

Also U2 = c'j,(/5fc, where the coefficients cJ, = {u2, fk) have at most polynomial 
growth ([25, Lemma 4.4]). Thus for every k G N there is mj. G N and > such 
that 

141 < 61.(1 + A,)'"'=. 

Define 

V = Vi ® U2 ^ ^ ® 'pk 

and let F be a cone in M"^ \{v ~ axis}. Then there is M > such that for all 
(A', A") e F we have that |A"| < A/|A'|. Therefore 

|c,c;,| < 6„,,(1 + A,)-"6U1 + AaO"^ 

< M"^6„j6'fc(l + Aj)-"+'"'= < C(M, j, fc)(l + A, + Afc)-"+'"^ 

and thus u is C°° on {Vti x X2) x F. 

Similarly, one can show that if xi G sing suppiti and U2 is smooth near X2, then 
u is C°° on (Xi X (sing supp(ii2))^) x L for all cones F in \ {a; — axis}. Therefore 

WF(ii) = (sing supp(iii) X sing supp(w2)) x 

U ((sing supp(?/i))'= x sing supp(M2)) x {y - axis} 
U (sing supp(ui) X (sing supp(u2))'^) x {a; — axis}. 

Example 8.6. We continue to assume that N ~2. Let x = {xi,X2) G X and let 
{(j, fc)} be a sequence such that the maps ipj ® ipk arc supported in a decreasing 
sequence of neighborhoods of x. Let {cjk} be any sequence of real numbers and 
define u = ^jk^j ® Vfc- Let f2 be an open subset of X such that there is an open 
neighborhood V of a; so that ilp]^ — ^- Then there are only a finite number of 
indices (j, k) with the property that supp (g) (^fc P| 17 7^ 0. It follows that u is C°° 
on X M^. Since was arbitrary, it follows that WF(u) C {x} x K^. Let F be a 
cone in R^. Define c^fe = 1 if {Xj^Xk) G F and set Cjk to be zero otherwise. Then 
for u = X Cjfc(pj ® (pk we have that WF(u) = {x} x F. 

More generally, let K be any compact subset of X^. Let {{j, fc)} be a sequence 
such that ipj <Si ipk is supported in a decreasing sequence of neighborhoods of K . 
Then if m = ^ '^jkfj (E)fk with {cjk} an arbitrary sequence, we have that WF(u) C 
K x R^. If F is a cone in R^ and we define Cjk as before, we obtain that WF(u) = 
K X Rl. 
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9. Pseudo-differential operators with variable coefficients 

The symbols of the operators studied in the previous section are independent 
of the X variable. These symbols are also known as constant coefficient symbols, 
and the corresponding operators are the constant coefficient operators. In this 
section we extend our study to operators whose symbols depend on both the A and 
X variables. We call them pseudo-differential operators with variable coefficients. 
The main difficulty in studying these operators is the fact that for a large class of 
fractals the domain of the Laplacian is not closed under multiplication [5]. This 
implies in our case that the symbolic calculus is not valid for the pseudo-differential 
operators with variable coefficients. Namely, the product of two symbols in the 
sense of Definition 9.1 is no longer smooth and it can not be a symbol of a pseudo- 
differential operator. Another consequence of this fact is that the kernel of these 
operators cannot be smooth. Nevertheless, we conjecture that these are Calderon- 
Zygmund operators in the sense of [31, 19]. We prove, in fact, that the kernels 
of these operators are continuous off diagonal and satisfy the correct decay off the 
diagonal. Moreover, we show that these operators are bounded on for all 

1 < q < oo. 

We can define the operators if X is either a compact fractafold without bound- 
ary based on a self-similar nested fractal or an infinite blow-up of K without 
boundary. However we prove the main properties in the case when X is compact. 
We continue to write V for the set of finite linear combinations of eigenvalues of 
—A with compact support and A for the spectrum of —A. We also assume that 
the heat kernel satisfies the upper estimate (2.2). 

Definition 9.1. For to G K we define the symbol class S"'" to consist of the smooth 
functions p : X x (0, oo) — C such that for each A; G N and j G N there is a positive 
constant Cj^k such that 



For A G A we set P\ to be the corresponding spectral projection. 

Definition 9.2. We define the operator class of pseudo-differential operators with 
variable coefficients 'i'DOm by 



for p G S"" and w G V. 

In the remainder of this section we will assume that Al is a compact fractafold 
without boundary. Recall that the set of eigenvalues is countable and the only 
accumulation point is oo. Let A = {Ai < A2 < . . . } be the set of eigenvalues of —A 
in increasing order and with repeated multiplicity. Let {(/'fc}feeN be an orthonormal 
basis of L^{n) consisting of compactly supported eigenfunctions. 

Theorem 9.3. Suppose that X is a compact fractafold with no boundary and p G 
S^. Then the operator Tu{x) = p{x,—/\)u{x) extends to a bounded operator on 



(9.1) 





AeA 



i2 
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Proof. Define 

fu{x,z) = Y^ J p{z,\)Px{x,y)u{y)dy. 



A6A 

Then 

\Tu{x)\ = \fu{x,x)\ < sup|rw(z,a;)| < ||(/ - A^)5(3TTTfM(-,a:)||2, 

z 

where the last inequality is the Sobolev inequality. So ||Tw||2 < — t^z^Tu\^. If 
p £ S'^ then AzP € L°° because X is compact. Then the result follows. □ 

Remark 9.4. The method of proof for the above theorem is adapted from [34]. 

Let K{x^ y) be the distributional kernel of T. Recall that we can not expect 
that K is smooth off the diagonal, because the domain of the Laplacian is not, in 
general, closed under multiplication. The best result one can expect is that K is 
continuous off diagonal. We prove that this is true in the following theorem. 

Theorem 9.5. Suppose that X is a compact fractafold with no boundary and p G 
S*^ . Then the operator Tu{x) ~ p{x, — A)ii(a;) is given by integration with a kernel 
that is continuous off the diagonal and satisfies the estimate 

(9.2) \Kix,y)\<Rix,y)-^. 

Proof. Since the map x t-^ p{x, A) is smooth for all A G R+, it follows that we can 
write 

oo 

(9.3) p{x,X)=J2mk{X)Mx), 

fc=0 

where mfc(A) = {p{-, X),(f)k), for all A G M+. We claim that fhk,n{X) := A^mfe(A) is 
a 0-symbol with constant coefficients, for all fc, n G N. For j > we have 

-(i^mfc.„(A) f ^^^n, / \ j / \ 

d)^ " J ^ QXIP^^-^^^kMx)d^iix) 

A^— p(x,A)(-A)"0fc(x)d/.(x) 
A^"(-A)"— p(a-,A)0fe(x)d/.(x), 
using the Gauss-Green formula. The estimates (9.1) imply that 



(9.4) 



^jd^mk,n{X) 



d\i 



< Q,n / \Mx)\dfl{x) < C,- „m(X)1/2. 



Thus TOfc^„(A) is a 0-symbol. Notice that the constants in (9.4) are independent of 
k. Now Theorem 4.5.4 of [21] (and the comments following it) implies that there is 
a positive constant c and a real number a such that 

(9.5) ||0fe||oo<cA^forallfc>O. 
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Let n > be such that the series S := J2k " converges. For u & V we can write 
then 



Tu{x) = ^^mk{Xj)(l)k{x)PxjU{x) 

k ^ 

where mfc^„(— A) is the pseudo-differential operator of order attached to rhk.n- 
Theorem 3.6 and estimates (9.4) imply that mfc_„(— A) is given by integration with 
respect to a kernel Kk^nix,y) that is continuous (and smooth) off diagonal and 
satisfies the estimate 

(9.6) \KkAx,y)\<C-Rix,y)-'', 
where C is a positive constant that is independent of k. Then 

Tu{x) = f ^^^^^Kk^n{x,y)u[y)diJ.{y). 
k ^fe 

Thus it suffices to prove that 



is well defined and continuous. Estimates (9.6) imply that 

\K{x, y)\<C- Rix, yy' Y^^^CS- i?(x, yy". 

k 

Thus K[x, y) is continuous off diagonal and (9.2) holds. □ 

Wc conjecture that the kernel K of a pseudo-differential operator p(x-, —A) of 
order is a Calderon-Zygmund operator in the sense of [31, Section 1.6.5]. Even 
though we are unable at this point to prove this claim, we can, nevertheless, show 
that these operators are bounded on L''{fi), for all 1 < g < oo. 

Theorem 9.6. Suppose that X is a compact fractafold with no boundary and p G 
5*^. Then the operator Tu{x) ~ p(.T, — A)w(x) extends to a bounded operator on 
L'^in) for alll < q < oo. 

Proof. Fix 1 < g < oo and let r such that 1/q+l/r = 1. Let u € V and let fhk^n as 
in the proof of Theorem 9.5. Recall that mfc,n is a symbol of order with constant 
coefficients, and 

Tu{x) = E ^T7r^^fc.n(-A)'u(a;), 



A: 



for all n > 0. Recall from inequality (9.5) that there is a such that Hi/ifelloo < cA^ 
for some positive constant c. Let n > such that the series Si — J2k •^l'" 
and 5*2 — J2k -^fc "^^^ ^^'^ convergent. Estimates (9.4) and Corollary 3.7 imply that 
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mk.n.{—^) is a bounded operator on L''(ii) with the bound independent on k. Then, 
Holder inequaUty imphes that 

q 

4>k{x) 



\Tu\\l 



< 



< 



E 



k 



dii{x) 



V fc 

E 



(t)k{x) mfe,„(-A)u(x) 



'/2 



dfi{x) 



\Mx)\ 



,nr/2 
k H 



E 



|mfe,„(-A)u(x)|9 



A 



nq/2 



E^ 



r(a— n/2) 



•^iEt;^ / l™/c,«(-A)u(a;)|''d/x(a;) 



< C{n)SiS2\\u\\l. 
Thus T extends to a bounded operator on L'^in). 



□ 
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